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The overlap hypercube fermion is a variant of a chirally symmetric 
lattice fermion, which is endowed with a higher level of locality than 
the standard overlap fermion. We apply this formulation in quenched 
QCD simulations with light quarks. In the p-regime we evaluate the 
masses of light pseudoscalar and vector mesons, as well as the pion decay 
constant and the renormalisation constant In the e-regime we present 
results for the leading Low Energy Constants of the chiral Lagrangian, 
E and F n . To this end, we perform fits to predictions by chiral Random 
Matrix Theory and by different versions of quenched Chiral Perturbation 
Theory, referring to distinct correlation functions. These results, along 
with an evaluation of the topological susceptibility, are also compared to 
the outcome based on the standard overlap operator. 



1 Overview 



QCD at low energy cannot be handled by perturbation theory Therefore 
it is often replaced by Chiral Perturbation Theory (%PT) as an effective 
theory pQ. deals with fields for the light mesons (occasionally also 

nucleons) which are the relevant degrees of freedom in that regime. Here 
we present non-perturbative results for QCD itself — the fundamental 
theory — at low energy, and we also establish a connection to the most 
important parameters in xPT. For such systems, chiral symmetry plays 
a central role, even though it is only realised approximately. 

In general, it is notoriously problematic to regularise quantum field 
theories with massless (or light) fermions in a way, which keeps track 
of the (approximate) chiral symmetry. For instance, in the framework 
of dimensional regularisation this issue is analysed carefully in Ref. [2]. 
However, in the framework of the lattice regularisation there was sub- 
stantial progress in this respect at the end of the last century (for a 
review, see Ref. |3]). At least for vector theories a satisfactory solution 
was found, which enables the simulation of QCD close to the chiral limit. 
But since this formulation is computationally demanding, production 
runs are limited to the quenched approximation up to now. 

In this work, we present simulation results employing a specific type 
of chiral lattice fermions, denoted as the overlap hypercube fermion 
(overlap-HF). In Section 2 we briefly review its construction and discuss 
some properties, which are superior compared to the standard overlap 
formulation. 

Section 3 is devoted to the regime, where the p-expansion of x?T jl] is 
applicable (p-regime). Here we present simulation results for bare quark 
masses in the range from 16.1 MeV to 161 MeV. We measure the masses 
of light pseudoscalar and vector mesons, and the quark mass according 
to the PCAC relation. The latter fixes the axial-current renormalisation 
constant Z4, which has a stable extrapolation to the chiral limit. We 
also present a matrix element evaluation of the pion decay constant F n , 
where, however, such an extrapolation is less stable. 

In a fixed volume V ~ (1.48 fm) 3 x (2.96 fm), we further decrease the 
bare quark mass to m q < 8 MeV, which takes us into the e-regime, i.e. the 
domain of the e-expansion [3] in x?T. The study of overlap-HFs in the 
e-regime, and also the comparison to the standard overlap fermions under 
identical conditions, is the main issue of this work, which we present in 
Section 4. In that regime, finite size effects and the topological sectors 
play an extraordinary role. 

Therefore, in Subsection 4.1 we first discuss the distribution of topo- 
logical charges — defined by the fermion index — and the corresponding 
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susceptibility. Then we address the Low Energy Constants (LECs) of 
the chiral Lagrangian, which parametrise the finite size effects. Hence 
they can be extracted even from the e-regime — i.e. from a relatively 
small volume — with their values in infinite volume, which are relevant 
in physics. In Subsection 4.2 we compare our results for the leading non- 
zero eigenvalue of the Dirac operator to the prediction by chiral Random 
Matrix Theory (RMT), which allows for a determination of the chiral 
condensate S. 

Finally we return to F n , which we evaluate in the e-regime by means 
of two methods. In Subsection 4.3 we fit the correlation function of 
the axial-currents to a prediction by quenched ;\RT, making use of the 
previously obtained values for and S. At last we work directly in 
the chiral limit (i.e. at zero quark mass) and consider the zero-mode 
contributions to the pseudoscalar correlation functions (Subsection 4.4). 
Again, a value for F n emerges from a fit of our data to a prediction 
by quenched ;\RT. However, the two versions of quenched that we 
apply to determine F n differ by a subtlety in the counting rules for the 
quenched terms in the e-expansion. 

We summarise and discuss our results in Section 5. Part of them 
were included previously in a Ph.D. thesis [H| and in several proceeding 
contributions [7j. 

2 The overlap hypercube fermion in QCD 

The lattice regularisation usually introduces a UV cutoff ir/a, where a 
is the lattice spacing. However, a block variable Renormalisation Group 
Transformation (RGT) generates a lattice action on a coarser lattice - 
say, with a spacing a' = 2a — which leaves the system unchanged, so 
that still the original cutoff matters jH]. Hence in this formulation the 
lattice artifacts are controlled by 27r/a', in contrast to the situation for a 
standard action. An infinite iteration of this blocking procedure leads in 
principle (for suitable RGT parameters) to a perfect lattice action, which 
is free of any cutoff artifacts. 

In practice, for most systems perfect actions can only be constructed 
in some approximation, such as a classical RGT step and a trun- 

cation of the long-range couplings is needed. Exceptions are, for instance, 
the quantum rotor JTTJ , the Gross- Neveu model at large N jT2] and the 
free fermion, where explicitly parametrised perfect actions are known 
[T3*l ITIj . Let us start from free Wilson-type fermions and iterate the 
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simple (blocking factor n) RGT 



Here x (a;') are the sites of the original, fine (blocked, coarse) lattice, 
populated by the spinor fields ip, ip (ip', ip'), with the lattice action S 
(5"), and x G x' are the sites x in the n d block with centre x' (in d 
dimensions, and n — 2,3 . . .). /i^Oisa real, dimensionless RGT 
parameter. After an infinite number of iterations, a perfect, free lattice 
action S*[ip,ip} emerges, with a Dirac operator D consisting of a vector 
term plus a scalar term, 

S*[7p, if)] = a d £ 4>xD xy il)y , D xy = 7 A1 p A1 (x - y) + \{x - y) , (2.2) 

x,y 

(where we denote the spacing on the final lattice again by a). For a finite 
RGT parameter /z, the scalar term A is non-zero. Then this action is local, 
since the couplings in p M and A decay exponentially at large distances 
\x — y\ even at zero fermion mass (x and y are lattice sites) ^Hj- This 
is satisfactory from the conceptual perspective, but in view of practical 
applications the couplings need to be truncated to a short range. A useful 
truncation scheme works as follows |15j : one first optimises the RGT so 
that the locality of the operator D is maximal, hence the truncation is 
minimally harmful. Then one constructs the perfect fermion action on 
a small, periodic lattice of size 3 4 (where we now specify d — 4). The 
couplings obtained in this way are finally used also in larger volumes. 
This is a truncation to a hypercube fermion (HF) operator D^f with 

supp[p^(x - y)) , supp[A(x - y)] C j x, y \x u - y v \ < a , W j . (2.3) 

The free HF still has excellent scaling properties [Tol I16j . 

At mass zero, the perfect fermions are chiral, since D solves the 
Ginsparg- Wilson relation (GWR) [T7j 

^75 + l5 D = -D l5 D (2.4) 
H 

for a simple RGT of the type (|2.1|) . with a transformation parameter 
The GWR is obvious for free perfect fermions, but it also persists 
under gauge interaction. Classically perfect Dirac operators solve the 



4 



GWR as well [TH], and the free HF operator D H f does so to very good 
approximation, as the spectrum shows [P3j . 

In our HF formulation of lattice QCD [13 1201 IE] 1 , the gauge field 
is treated rather ad hoc, since it is the fermionic part which is most 
delicate. Hence we use the standard Wilson gauge action, and the Dirac 
operator is gauged as follows: if the lattice spinors ip x and ip y are coupled 
by Dnp,xy — i- e - if they are located in the same lattice unit hypercube - 
we connect them by the shortest lattice paths and multiply the compact 
link variables on these paths. The normalised sum of these link products 
is denoted as a hyperlink between ip x and ip y ( a fully explicit description 
is given in Refs. |2TH IB]). 

In addition, we replace the simple links in the above procedure by fat 
links, Ufj,(x) — > (1 — at)Un(x) + | Yl [staples] , where a is a parameter to 
be optimised. 

This HF is nearly rotation invariant, but its mass is subject to a 
significant additive renormalisation rj^j . Criticality can be approximated 
again by amplifying all the link variables by a suitable factor u > 1. At 
last, in the vector term only they are also multiplied by an extra factor 
v < 1; this hardly alters the critical point, but it improves chirality, i.e. 
it moves L>hf closer to a GWR solution. The fat link, e.g. with a in the 
range 0.3 . . . 0.6, helps further in this respect [22 El- 

An optimisation of that kind was performed previously for quenched 
QCD at (3 = 6 [23 J. Here we present the more difficult construction at 
f3 = 5.85, which corresponds to a lattice spacing of a ~ 0.123 fm. 2 As a 
compromise between various optimisation criteria, we chose the following 
parameters [HI 

u = 1.28, w = 0.96, a = 0.52, // = 1 , (2.5) 

and the HF couplings are still those identified for the massless, truncated 
perfect, free fermion ^3]. The physical part of the spectrum of this .Dhf 
for a typical configuration at — 5.85 is shown in Fig. ^ We also mark 
some of the eigenvalues with maximal real part; they reveal a striking 
difference from Wilson fermion spectra, which extend to much larger real 
parts. Moreover, we include the low eigenvalues for an alternative set of 
parameters, which moves the HF close to criticality. In that case we 
recognise a good approximation to the circle in the complex plane with 
centre and radius /i = 1, which characterises an exact solution to the 
GWR (|2.4|) . However, for the practical purposes in this paper, we prefer 
the parameters in eq. (J2.5)) (the advantages are a faster convergence and 

^or recent applications in thermodynamic systems, see Ref. |21|. 
2 Throughout this work, we refer to the Sommer scale [2, for physical units. We 
do not keep track of possible errors in that scale. 
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a better locality if the HF operator is inserted in the overlap formula, 
which we are going to discuss next). 




0.5 1 1.5 2 

a Re X 

Figure 1: The low lying part of the spectrum o/D ov hf ; plus eigenvalues 
with maximal real part, for a typical configuration at (3 = 5.85 on a 8 4 
lattice. We show one set of parameters which provides a good approxima- 
tions to criticality and to the spectral Ginsparg- Wilson circle. However, 
for practical purposes we prefer the parameterisation \2.5)) . which leads 
to the other set of eigenvalues shown in this Figure. 

Of course, due to the truncation and the imperfect gauging, the scal- 
ing and chirality of this HF are somewhat distorted. At least chirality 
can be corrected again by inserting .Dhf into the overlap formula [23] 

D av =?-(l + A/VAtA) , A:=D -^. (2.6) 

This formula yields a solution to the GWR ()2.4j) . provided that D is 
a 75-Hermitian lattice Dirac operator, D\ = 75/^075- Another condition 
is that a value for /i must be available, which separates — for typical 
configurations at the gauge coupling under consideration — the (nearly) 
real part of the spectrum of Do in a sensible way into small (physical) 
eigenvalues and large eigenvalues (at the cutoff scale). Thus the physical 
eigenvalues are separated from the doublers, and we obtain the correct 
number of flavours. 75-Hermiticity holds for example for the Wilson 
operator _D W and for -Dhf, and for instance at /3 = 6 good values for /i 
can easily be found in both cases • 

D = D w is the kernel of the standard overlap operator formulated by 
H. Neuberger |25j : we denote it by D-^. In this case, -Dw changes drasti- 
cally as it is inserted into the overlap formula (|2.(jjl . In contrast, D^f is 
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already an approximate Ginsparg- Wilson operator, hence the transition 
to the corresponding overlap operator, _D H f — ¥ ^ovUF, is only a modest 
modification, 

D ovUF w D UF , (2.7) 

particularly in view of the eigenvalues with large real parts. In both 
cases, a lattice modified but exact chirality is guaranteed |2Zj through 
the GWR. The correct axial anomaly is reproduced in all topological 
sectors for the standard overlap fermion [28., as well as the overlap-HF 
[29J. But other important properties are strongly improved for -D ov hf 
compared to D^, due to the construction of Duf and relation ([2.7)1 . 3 

Let us first consider the condition number of the Hermitian operator 
A^A, i.e. the argument of the square root in the overlap formula (|2.6|) . 
In our simulations, we project out the lowest 30 modes of A^A — which 
are treated separately — and we approximate the remaining part by a 
Chebyshev polynomial to an absolute accuracy of 10~ 12 . The polynomial 
degree which is needed for a fixed precision is proportional to the square 
root of the condition number 

C31 := max , with A max (A31) : maximal (31 st ) eigenvalue of A^A. 

(2.8) 

The values for (C31), given in Table 1, show that the required polyno- 
mial degree for -D ov hf is a factor 5 lower. The computational effort 
is roughly proportional to this degree (the computation of the lowest 
modes of A^A and their separate treatment are minor issues in terms of 
the computing time). On the other hand, the use of the HF kernel is 
computationally more expensive than Dyj by about a factor 15, so that 
an overhead by a factor « 3 remains. We hope for this factor to be more 
than compensated by the virtues of -D ov hf, in particular by an improved 
scaling due to the perfect action background of Dhf and relation ()2.7|) . 
That property is very well confirmed for the free fermion JT5| and for 
the 2-flavour Schwinger model (with quenched configurations, but mea- 
surement data re- weighted by the fermion determinant) [30]. But the 
corresponding scaling behaviour in QCD is not explored yet. 

An important aspect that we do explore here is the locality. For this 
purpose, we put a unit source at the origin, r) x = 8 X) o, and measure the 
expectation value of the function j2E] 

/max(ri) := m , ax [\\D OVtyx (U) Vx \\ I \\y\U = n} . (2.9) 



3 The concept of improving properties of D ov beyond chirality by choosing a more 
favourable kernel than Dyv was suggested in Ref. 19 and meanwhile implemented in 



several variants [23 E2 EH HI3 E2 • 
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P 




-DovHF 


(Cll) 


5.85 


4266(194) 


179(10) 


(on) 


5.85 


1723(46) 


73(2) 


<<*i> 


5.85 


1149(15) 


49.2(8) 


(7loc) 


5.85 


0.63(1) 


0.73(1) 


(7loc) 


5.7 


0.56(2) 


0.76(5) 


(7loc) 


5.6 




0.53(3) 



Table 1: The characteristic indicators for the kernel condition number 
and for the locality of the overlap operators (standard) and -D ov hf 
(described in the text), on a 12 3 x 24 lattice. 

Following the usual convention we take here the distance in the taxi 
driver metrics ||y||i := ^ A comparison, still at (3 = 5.85, is given 

in Fig. El (on top) and Table [TJ we see a clear gain for -D ov hf, i-e. a decay 
of fiji) oc exp(— 7ioc^i) (beyond short r%) where 7i oc is larger for D oviiF 
than for D^. If we refer to the Euclidean distance r = \y\, the maximum 
is taken over different sets, but the resulting decay is still exponential. 
In the lower plot of Fig. |21 we see that the decay for -D ov hf is not only 
faster, but it also follows in a smoother way the exponential shape in r, 
which hints at a good approximation to rotation symmetry. 

Locality is a vital requirement for a safe continuum limit in lattice 
field theory This poses a limitation on the overlap operator (J2.6|) . 
since for decreasing /3, at some point an exponential decay is not visible 
any longer. 4 

Fig- El shows the situation at (3 — 5.7 (i.e. a ~ 0.17 fm), where both 
operators are still local (with optimal parameters \x and u for resp. 
-Dovhf), but again the locality of -D ov hf is superior. As we push further 
to (3 — 5.6, we do not observe an exponential decay for anymore, but 
-Dovhf is still manifestly local (even without adjusting the parameters 
in eq. (|2.5|h except for u). Therefore, the overlap-HF provides chiral 
fermions on coarser lattices. 

3 Applications in the ^-regime 

To handle strong interactions at low energy, %PT replaces QCD by an 
effective Lagrangian. It involves the terms allowed by the symmetries, 

4 The question, in which sense locality still persists in a conceptual sense is discussed 
in Ref. jSHJ. In any case, if no exponential decay is clearly visible, the operator should 
certainly not be applied in simulations. 
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Figure 2: TTie locality of different overlap fermions, measured by the 
maximal impact of a unit source rj x over a taxi driver distance r\ ( on 
top), and over a Euclidean distance r (below), at [3 = 5.85. We see that 
the overlap-HF is clearly more local. The upper plot also shows that the 
slope hardly changes if we proceed from a 10 4 to a 12 3 x 24 lattice. (The 
increasing slope beyond t\ = 18 is a consequence of the anisotropy of the 
12 3 x 24 lattice.) The lower plot illustrates additionally that the overlap- 
HF decay follows closely an exponential in the Euclidean distance, which 
indicates a good approximation to rotation symmetry. 



which are ordered according to a low energy hierarchy. To its leading 
order, the effective Lagrangian of takes the form 



C eS [U] = ^Tr\d IM U(xYd tl U(x) 



Tr 



U{x) + U( X y +... . (3.1) 
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Figure 3: The locality of different overlap fermions, measured by the 
maximal impact of a unit source rj x over a taxi driver distance r\, on a 
12 3 x 24 lattice. For the overlap-HF we see a superior locality at (3 = 5.7. 
Its locality persist up to P = 5.6, where no exponential decay is visible for 
the standard overlap formulation. 



Throughout this work, we assume a degenerate quark mass m q for all 
flavours involved. The LECs are the coefficients attached to the terms 
in this expansion; in the leading order F n and E. The effective theory 
does have its predictive power, but the values of the LECs cannot be 
determined within %PT. To this end, one has to return to QCD as the 
fundamental theory, and it is our main goal in this work to investigate 
how far the LECs can be obtained from lattice QCD simulations. 

In xPT, one usually considers a finite spatial box, which has a volume 
L 3 , and one expects the meson momenta p to be small, so that 

p ~ — < 4nF n . (3.2) 

Note that the term AttF^ takes a role analogous to Aqcd- 

%PT then deals with a perturbative scheme for the momenta and 
masses of the light mesons, with appropriate counting rules. The most 
wide-spread variant of %PT assumes the volume to be large, I > ( = 
I/run (where £ is the correlation length, given by the inverse pion mass). 
Then the p- expansion |3] can be applied. This is an expansion in the fol- 
lowing dimensionless ratios, which are expected to be small and counted 
in the same order, 

TV \ \ ' (3,3 ^ 

J^^tt aqcd ^qcd 
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In this Section we present our simulation results in the p-regime. We 
apply the overlap-HF described in Section 2, at (3 = 5.85 on a lattice of 
size 12 3 x 24, which corresponds to a physical volume of V ~ (1.48 fm) 3 x 
2.96 fm. We evaluated 100 propagators for each of the bare quark masses 



(in physical units: 16.1 MeV . . . 161 MeV), using a Multiple Mass Solver. 
We will see that the smallest mass in this set is at the edge of the p-regime 
- even smaller quark masses will be considered in Section 4. 
We include m q in the overlap operator ()2.6|) in the usual way, 



which leaves the largest real overlap Dirac eigenvalue invariant. m q rep- 
resents the bare mass for the quark flavours u and d. 

We first evaluate the pion mass in three different ways: 

• m n,pp is obtained from the decay of the pseudoscalar correlation 
function (P(x)P(O)), with P(x) = ipxjsipx- This is the most ob- 
vious method, but it is not the best one in this case, as we will 
see. 

• m Tr,AA is extracted from the decay of the axial-vector correlation 
function (At(x)A 4 (0)), with A±{x) = ^xlhl^x- 

• fn-K,pp-ss is obtained from the decay of the difference (P(x)P(O) — 
S(x)S(0)), where S(x) = i[) x ip x is the scalar density. This sub- 
traction is useful especially at small m q , where configurations with 
zero-modes ought to be strongly suppressed by the fermion deter- 
minant. In our quenched study, this suppression does not occur as 
it should, but the above subtraction in the observable eliminates 
the zero-mode contributions, which are mostly unphysical. 

We present the results in Fig.HJand Tabled (the latter collects all the 
results of this Section). The pion masses follow to a good approximation 
the expected behaviour oc m q . Deviations show up at the smallest 
masses, where we observe the hierarchy 



in agreement with Ref. [TO]- This shows that the scalar density subtrac- 
tion is in fact profitable, since it suppresses the distortion of the linear 
behaviour down to the lightest pion mass in Fig. HI 



am q = 0.01, 0.02, 0.04, 0.06, 0.08 and 0.1 




(3.4) 



m^pp > m nAA > m 7ri pp_ S5 , 



(3.5) 



m n ,p P -ss(am q = 0.01) = (289 ± 32) MeV . 



(3.6) 
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That mass corresponds to a ratio L/£ ~ 2, hence around this point we 
are indeed leaving the p-regime. Based on the moderate quark masses in 
Fig. HJ we find an impressively small intercept in the chiral extrapolation, 

m^p P -ss{m q -> 0) = (-2 ± 24) MeV . (3.7) 

On the other hand, at our larger m q values the hierarchy changes due to 
the interference of the scalar correlator. 




am. 



Figure 4: The pion mass evaluated from overlap HFs in the p-regime in 
three different ways, as described in Section 3. 



Due to quenching, one expects at small quark masses a logarithmic 
behaviour of the form 

o 

Qfffl 

= C\ + C 2 In am„ + C 3 am q , (Ci, C 2 , C 3 : constants) . (3.8) 

rn q 

Corresponding results are given for instance in Refs. [HU Fig. 13.81 
shows the fits of our data to eq. (|3.8p . m n ^A is best compatible with this 
rule (this property was also hinted at in Ref. [H3]). At least the deviation 
for m^^pp^ss (although inside the errors) could be expected, since the 
scalar subtraction alleviates the quenching artifacts, which give rise to 
the logarithmic corrections according to formula ([3.8)1 . 

In Fig. |U1 (and Table EJ) we consider the vector meson mass m p based 
on the same 100 configurations. A chiral extrapolation leads to 

m p = (1017 ± 40) MeV . (3.9) 
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0.1 


0.08 


0.06 


0.04 


0.02 


0.01 




0.491(5) 


0.438(5) 


0.380(5) 


0.314(4) 


0.248(5) 


0.210(11) 




0.483(4) 


0.432(5) 


0.378(5) 


0.317(6) 


0.240(7) 


0.189(9) 


am 7T,pp-ss 


0.513(7) 


0.453(8) 


0.381(9) 


0.313(11) 


0.221(15) 


0.181(20) 


am p 


0.745(7) 


0.717(8) 


0.694(10) 


0.677(15) 


0.674(32) 


0.672(55) 


ampcAc 


0.089(3) 


0.070(3) 


0.052(3) 


0.035(3) 


0.017(2) 


0.0087(13) 


Z A 


1.13(4) 


1.14(5) 


1.15(6) 


1.16(9) 


1.16(14) 


1.14(18) 




0.117(2) 


0.114(3) 


0.110(3) 


0.105(4) 


0.103(4) 


0.098(4) 


a F n) pp-ss 


0.120(2) 


0.114(3) 


0.109(3) 


0.104(5) 


0.010(10) 


0.082(12) 



Table 2: Our results in the p-regime for the pion mass, the p-meson mass, 
the PC AC quark mass, the renormalisation constant and the pion 
decay constant F^, as described in the text and plotted in the Figures of 
Section 3 (Figs. . . 0). These results are obtained from 100 propagators 
at each of the bare quark masses m q . All the dimensional numbers in this 
Table are given in lattice units at a ~ 0.123 fm = (1610 MeV) -1 . 
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Figure 5: The pion masses of Fig. ^| fitted against the form \8. which 
is expected for the logarithmic quenching artifacts in the absence of an 
additive mass renormalisation. 



This agrees well with a study using the standard overlap operator on 
the same lattice ^jg^. It is quite large, however, not only in view of 
phenomenology, but also compared to other quenched results in the lit- 
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erature, which were summarised in Ref. |3Z|. 5 

On the other hand, if we insert our measured ratio m pseudoscalar /m vector 
at m q = 0.01 into the phenomenologically motivated interpolation for- 



mula presented in Ref. [38], then we arrive at a significantly lower value 
of m p » 789 MeV. 

Fig. [7] shows the quark mass obtained from the PCAC relation, 



?™PCAC 



Eg«Mto)g(o)) 
E,(Pt(x)P(o)) 



(3.10) 



which follows closely the bare quark mass m q (we use a symmetric nearest- 
neighbour difference for 84). As a consequence, the axial-current renor- 
malisation constant 

(3.11) 



"^PCAC 

is close to 1, see Fig.|H]and Table^J A chiral extrapolation leads smoothly 
to 

Z A = 1.17(2) , (3.12) 

which is in contrast to the unpleasantly large values found for the stan- 
dard overlap operator: e.g. at the same (3 = 5.85 and ji = 1.6 (a preferred 

5 Note that the difference from the data by other groups — obtained with various 
(mostly non-chiral) lattice formulations — is approximately constant in m T over the 
range shown here, hence finite size effects can hardly be blamed for it. 
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Figure 7: The PC AC quark mass, evaluated for -D ov hf according to eq. 
\3.1U\) . from stable plateaux in the time direction. We find mpcAC values 
close to m q , in contrast to the results with the standard overlap operator. 
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Figure 8: The axial-current renormalisation constant determined 
from eq. H3.ll)) . which is found to be close to 1 for the overlap-HF. Za 
is plotted against PP _ S s- 



15 



value in that case) it amounts to Z& ~ 1.45 (UHl EH1 HO] , and (somewhat 
surprisingly) at (3 = 6, /i = 1.4 it even rises to ~ 1.55 [30] ■ Accord- 
ing to Ref. [22] 5 the fat link could be especially helpful for the property 
Za ~ 1, which is favourable for a connection to perturbation theory. 

As a last observable in the p-regime, we measured the pion decay 
constant by means of the (indirect) relation 

F, = ^f |(0|P|7r)| , (3.13) 

based on P(x)P(0), or based on P(x)P{0) - S(x)S(0) (in eq. lETTSJl this 
affects both, the denominator and the pion state). 6 The results are 
given in Fig. Eland Table El In particular the value at am q = 0.01 (the 
lightest quark mass in this plot) is significantly lower for the case of 
the scalar subtraction. Hence this pushes the result towards the chirally 
extrapolated phenomenological value of ~ 86 MeV 




0.35 



Figure 9: The pion decay constant based on a matrix element evaluation 
in the p-regime — given by eq. — using the overlap-HF. 



However, a chiral extrapolation based on these data for F w would be 
risky — apparently we are too far from the chiral limit for this purpose. 
An extrapolated value of F n) pp would come out clearly too large, as it 
is also the case for [3^]. But in particular the instability of F n ,pp-ss 



6 A recent alternative attempt to evaluate LECs from simulations in the p-regime 
was presented in Ref. [^Tj . 
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at am q = 0.04, 0.02, 0.01 calls for a clarification by yet smaller quark 
masses. We did consider still much smaller values of m q in the same 
volume. As the results for the pion masses suggest, we are thus leaving 
the p-regime. For the tiny masses am q < 0.005 we enter in fact the e- 
regime, where observables like F n have to be evaluated in a completely 
different manner. This is the subject of Section 4. 



4 Applications in the e-regime 

The e-regime of QCD is characterised by a relatively small volume, i.e. 
the correlation length £ exceeds the linear box size L. On the other hand, 
the box still has to be large compared to the scale given by 4:irF n , as in 
eq. (|3.2j) . so that the partition function is saturated by the lowest states 
only — higher states (above about 1 GeV) do not contribute significantly. 
This amounts to the condition 

— > L > — — . 4.1 

In such a box, the p-expansion of %PT fails, in particular due to the 
dominant role of the zero-modes. However, the latter can be treated 
separately by means of collective variables, and the higher modes - 
along with the pion mass — are then captured by the e-expansion [5]. 
One now counts the ratios 

p2 1 (4.2) 



Aqcd A^ cd (LF w ) 2 

as small quantities in the same order. 

The Haar measure for the fields U(x) G SU (Nf) in the coset space of 
the chiral symmetry breaking, 

SU(N f ) L ® SU(N f ) R -> SU(N f ) L+R (for N f flavours) , (4.3) 

is worked out in Ref. [13]. The corresponding non-linear O(N) a- model 
(with a symmetry breaking to 0(N — 1)) in a small box was studied with 
the Faddeev- Popov method to two loops jH], and with the Polyakov 
functional measure to three loops 

This setting cannot be considered a physical situation. Nevertheless 
there is a strong motivation for its numerical study: the point is that 
the finite size effects are parametrised by the LECs of the effective chiral 
Lagrangian as they occur in infinite volume, hence the physical values of 
the LECs can (in principle) be evaluated even in an unphysically small 
box. 
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We recall that the effective chiral Lagrangian includes all terms which 
are compatible with the symmetries, ordered according to suitable low 
energy counting rules, in this case the counting rules (|4.2j) of the e- 
expansion. We further presented in Section 3 the LECs as coefficients 
of these terms, for instance F n and £ in the leading order. Their de- 
termination from QCD — the fundamental theory — at low energy is 
a challenge for lattice simulations. As a test, such a lattice determina- 
tion in the e-regime has been performed successfully for the 0(4) cr-model 
(which describes chiral symmetry breaking with two flavours) some years 
ago j3H], but in QCD this method jUj had to await the advent of chi- 
ral lattice fermions. Unfortunately, the quenched results for the LECs 
are affected by (mostly logarithmic) finite size effects |48j . so that the 
final results by this method still have to wait for the feasibility of QCD 
simulations with dynamical, chiral quarks. 

A peculiarity of the e-regime is that the topological sector plays an 
essential role jlH] : if one measures observables in a specific sector, the ex- 
pectation values often depend significantly on the (absolute value of the) 
topological charge in this sector. In particular for the evaluation of the 
LECs, a numerical measurement inside a specific sector and a confronta- 
tion with the analytical predictions in this sector is in principle sufficient. 
This requires the collection of a large number of configurations in a spe- 
cific sector. The "topology conserving gauge actions" (5U] are designed 
to facilitate this task. However, here we stay with the Wilson gauge ac- 
tion, which allows us to investigate also the statistical distribution of the 
topological charges, which we are going to address next. 

4.1 The distribution of topological charges 

A priori, it is not obvious how to introduce topological sectors in the 
set of lattice gauge configurations. However, if one deals with Ginsparg- 
Wilson fermions, a sound definition is given by adapting the Atiyah- 
Singer Theorem from the continuum and defining the topological charge 
by the fermionic index v |18j . 

topological charge = v := n + — n_ , (4-4) 

where n± is the number of zero-modes with positive/negative chirality. 
We remark that these numbers are unambiguously determined once a 
Ginsparg- Wilson Dirac operator is fixed (and that in practice only chi- 
rality positive or chirality negative zero-modes occur in one configura- 
tion 7 ). However, for a given gauge configuration, the index for different 

7 This property even holds in cases where "cooling" deforms the configuration into 
a form, where a semi-classical picture suggests the presence of topological objects 
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Ginsparg- Wilson operators does not need to agree. Albeit the level of 
agreement should be high for smooth configurations, i.e. it should — and 
it does 8 — increase for rising values of f3. 

In our study, still at j3 — 5.85 on a 12 3 x 24 lattice, we compared the 
charges for the overlap-HF operator described in Section 2 and for the 
standard overlap operator at fi = 1.6. As an example, the histories 
of about 200 indices for the same configurations are compared in Fig. 
ITU1 Of course, these two types of indices are considerably correlated, but 
only 41% really coincide. We obtained a mean deviation of 

(l^ovHF — vn\) = 0.80(2) , (4.5) 

and we observed over more than 1000 configurations a maximal index 
difference of |z/ ov hf — ^n| =5. Still, the similarity is of course much closer 
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X 
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configuration number 

Figure 10: Index histories for D oviiF (see Section 2) and for _D N (at 
fi = 1.6) for the same set of configurations. 

than the accidental agreement for independent indices, since they follow 
essentially the expected Gaussian distributions, with a width m 3.3, see 
Fig. [TT] and Table El This width fixes the topological susceptibility 

with opposite chiralities Regarding the fermion index, a cancellation happens 

for instance for free Ginsparg- Wilson fermions, but in a realistic gauge background 
such an unstable constellation is very unlikely. 

8 For instance, we observed at j3 = 6.15 on a 16 4 lattice that the index of Z?n is 
very stable as /i rises from 1.3 to 1.7; this changes less than 2% of the indices. 
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Figure 11: The histogram of the overlap-HF indices (on the left) and for 
the standard overlap indices. In both cases, 1013 indices are included. 



Xtop = ^ 2 >, (4.6) 

which is of importance to explain the heavy mass of the rj' meson |52) . A 
discussion of that point, as well as a measurement of Xto P based on 
indices on L 4 lattices with a continuum extrapolation, is given in Refs. 

In Fig. H21 and Table EH we present our results with -D ov hf and D N 
on the lattice used so far, plus a result for at (3 = 6 in the same 
physical volume (lattice size 16 3 x 32). We also mark the continuum 
extrapolation according to Ref. [53], which is fully consistent with our 
results. This value for Xtop is compatible with the Witten-Veneziano 
scenario that much of the r/ mass is generated by a U(l) anomaly. 

In principle, the charge histogram could give insight into the possi- 
bility of a spontaneously broken parity symmetry in QCD, which is not 
fully ruled out 54]. This question was also studied in Ref. [55] with a 
different definition of topological charges on the lattice. Here we observe 
in all cases that the number of neutral configurations is about half of the 
corresponding number with \v\ — 1 (see Table Based on this observa- 
tion, it cannot be decided if the charge distribution (at small favours 
a precise Gaussian, or a double peak structure (or something else). Hence 
the fate of parity symmetry remains open. 



4.2 Determination of the chiral condensate £ 

Chiral RMT conjectures predictions for the low lying eigenvalues, or- 
dered as A n , n = 1,2,... (excluding possible zero eigenvalues) of the 
Dirac operator in the e-regime (for a review, see Ref. [5B]). More pre- 
cisely, the conjectured densities are functions of the dimensionless vari- 
ables TV\ n) where E is the chiral condensate in the effective Lagrangian 
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total ^ of confs. 


lOlo 


lOlo 


c n^ 
506 




in 01 i r\ a^7 


11 /in i n r 1 
11. 4y ± U.Ol 


in /i o _i_ n 
lU.4y ± U.DD 


4 

Xtop r 


0.071(3) 


0.076(3) 


0.069(4) 


^ of confs. with v = 


118 


124 


59 


# of confs. with z/ = 1 


238 


237 


115 


# of confs. with z/ = 2 


210 


192 


95 


^1/3 f rom sec tor = 


298(4) MeV 


301(4) MeV 


279(7) MeV 


a(Ai)|„=o 


0.0069(4) 


0.0067(4) 


0.0059(4) 


a(Ai)|„=i 


0.0130(4) 


0.0136(4) 


0.0093(3) 


a(Ai)|^=2 


0.0184(5) 


0.0193(5) 


0.0130(5) 


a(Ai)| I /= 3 


0.0247(6) 


0.0242(8) 


0.0155(6) 


a(\ 1 )\ u =4 


0.0293(9) 


0.0312(10) 


0.0215(9) 


a(Ai)|^= 5 


0.0338(12) 


0.0360(13) 


0.0246(17) 


S 1 /^ from (Ai)| y | = o...5 


290(6) MeV 



Table 3: Our results for the topological susceptibility and for the chiral 
condensate in the e-regime for a fixed physical volume V = (1.48 fm) 3 x 
2.96 fm. We consider two types of overlap operators (D miiF as described 
in Section 2, and the standard overlap operator at fi = 1.6), and 
two lattice spacings. We first give the total statistics and the resulting 
topological susceptibility (see Subsection J^.l) in a dimensionless form, 
with ro = 0.5 fm (according to the Sommer scale WM)- 
Below we give separately the statistics in the sectors \v\ = 0, 1 and 2. 
We extract the chiral condensate £ from the density of the lowest Dirac 
eigenvalue in the neutral charge sector, which is most reliable since it 
involves the smallest values of z, see Subsection 4-2 and Fias. f73l \T\ 



As an alternative we considered the mean values of the lowest non-zero 
Dirac eigenvalues (Ai) in the sectors \v\ = ... 5. For the value of £ in 
the last line all these results for (Ai)^ match the RMT predictions, see 
Fig.M 
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Figure 12: The topological susceptibility measured by indices of -D ov hf 
and of _Dn; ^ n & volume V = (1.48 fm) 3 x 2.96 fm ; with two different 
lattice spacings a. Our data agree well with the continuum extrapolation 
reported in Ref. 

(13. lj) . Here we focus on the variable z := EVA^p, where A^p emerges 
from the leading non-zero eigenvalue Ai if the spectral circle of the over- 
lap operator is mapped stereographically on the imaginary axis, Ai,p = 
|Ai/(l-aAi/2/x)|. 9 

These RMT predictions depend on \u\, the absolute value of the topo- 
logical charge. Here we make use of the explicit formulae [HZ] for the den- 
sity of the first non-zero (re-scaled) eigenvalues in the sectors \v\, which 
we denote by p[ (z). For the lowest eigenvalues, the particular den- 
sity Pi was first confirmed by staggered fermion simulations (results are 
summarised in Ref. jSH]), but the charged sectors yielded the very same 
density, in contradiction to RMT. The distinction between the topologi- 
cal sectors was first observed to hold for to a good precision [HE], if 
the linear box size exceeds a lower limit of about L> 1.1 fm (of course, 
the exact limit depends on the criterion). 10 Once the predicted density 
Pi"^ is well reproduced, we can read off the value of S, which is the only 
free fitting parameter for all topological sectors. It is most instructive 

9 Alternatively, one could simply consider |Ai|, which are eigenvalues of 75-D, but 
for the small eigenvalues that we deal with the difference is not of importance. 

10 Meanwhile, a topological splitting was also observed to set in for staggered fermion 
if the link variables are strongly smeared |59j. 
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to plot the cumulative densities JOU], which we show in Figs. Upland ITU 
We compare here the predictions to the eigenvalues X^p, which we mea- 
sured in various topological sectors. The statistics involved in each case 
is included in Table EJ This Table also displays the S values obtained 
in the sector v = 0, which we consider most reliable, since it deals with 
the lowest eigenvalues resp. energies. As a theoretical bound, one often 
refers to the Thouless energy F^/(T,yV), below which these predictions 
should hold. In our case, it translates into ZThouiess^l, but the eigen- 
value distributions follow the chiral RMT behaviour up to larger z values. 
Clearly, in this range it is the neutral sector {y = 0) which contributes in 
a dominant way, but Fig. 1141 shows that in the case of (for both values 
of j3) the charged sectors \u\ = 1 and 2 alone would favour a different S 
value. This ambiguity also occurs for smeared staggered fermions jSH]- 
In the case of the -D ov hf ; however, a unique S works well for all the three 
sectors \u\ — 0, 1, 2, up to about z ~ 3, as we see from Fig. [TBI 



12 3 x24, (3=5.85 




z 



Figure 13: The cumulative density of the lowest Dirac eigenvalue X± t p of 
the overlap-HF operator, in the topological sectors \u\ =0,1 and 2. We 
compare the chiral RMT predictions to our data for z = EVA^p with 
X 1 / 3 = 298 MeV — the optimal value in the neutral sector (v = 0). This 
value works well up to z w 3 for all topological sectors. 

As an alternative approach to test the agreement of our data with the 
chiral RMT, and to extract a value for S, we now focus on the mean val- 
ues of the leading non-zero Dirac eigenvalues Ai in all the charge sectors 
up to \v\ = 5. In physical units, the results (Xi t p) agree remarkably well 
for the different overlap operators and lattice spacings — see Fig. [Tol 
although this consideration extends beyond very low energy. Each single 
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Figure 14: The cumulative density the lowest Dirac eigenvalue X^p of the 
standard overlap operator at (3 = 5.85 (on top) and (3 — 6 (below), 
in the topological sectors \u\ = 0, 1 and 2. We compare the chiral RMT 
predictions to our data for z = T.VX^p with S 1 / 3 = 301 MeV (on top), 
and with S 1 / 3 = 279 MeV (below) — the optimal values in the neutral 
sector (v = 0). Considering also \u\ = 1,2 would decrease (increase) S 
at (3 = 5.85 ((3 = 6) which is the trend towards the result with the method 
of Fig. [H 

result for (Ai p)^ can then be matched to the RMT value for a specific 
choice of S. Amazingly, all these 18 results are in agreement with the 
RMT if we fix 

E = (290(6) MeV) 3 , (4.7) 
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as Fig. El also shows. This value is between the results obtained from 
the eigenvalue densities at v = alone, and we recognise from Figs. El 
and El that this is the trend if we take the charged sectors into account. 

A renormalisation procedure for £ obtained in this way is discussed 
in Ref. [HIJ- However, we will only use this quenched lattice result as a 
fitting input in Section 4.3, so here we stay with the bare condensate E 
for our fixed lattice parameters. 



60 
50 


i 


i i i 

overlap-HF, (3=5.85 x 
standard overlap, (3=5.85 i — i — i 
standard overlap, (3=6 s b 
RMT at £ = 290(6) MeV • 


i 

1 
• 


i 

I" 

"fP_ 


40 




* 






30 




1 






20 




* 






10 


* 
i 


i i i 


i 


i 







1 2 3 


4 


5 



|v| 



Figure 15: The mean values of the first non-zero Dirac eigenvalue in the 
charge sectors \v\ = . . . 5, in physical unit. All the measured results 
agree with chiral RMT if we choose S 1 / 3 = 290(6) MeV. 



4.3 Evaluation of F n based on the axial-vector cur- 
rent correlator 

As we mentioned before, QCD simulations with chiral quarks can only be 
performed in the quenched approximation for the time being. In order to 
relate simulation results to the effective low energy theory, we therefore 
refer to quenched xPT. In that framework, mesonic correlation functions 
were evaluated to the first order in Refs. [021 ESj • It turned out that the 
vector current correlation function vanishes; this property actually ex- 
tends to all orders [64]. The scalar and pseudoscalar correlators involve 
already in the first order additional, quenching specific LECs, which ob- 
struct the access to the physical LECs in the Lagrangian 1)3. lj) . Therefore 
we first focus on the axial-vector correlator, which only depends on E and 
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Fn in the first order (as in the dynamical case). In particular we are go- 
ing to compare our data to the quenched ;^PT prediction in a volume 

V = L 3 x T 163 , 



Z\.(A A {t)AM) 





V 




E z q I u (z q )K u (z q ) + !„+! + 




(4.8) 



where 



A 4 (t) = a 3 J~]i(){t, f)7574^(t, x) 



(f >0) 



(4.9) 



X 



is the bare axial- vector current at 3- momentum p = 0. This formula 
applies to the topological sectors of charge ±u. I v and K v are modi- 
fied Bessel functions, and z q := HVm q (in analogy to the variable z in 
Subsection 4.2). 

It is remarkable that this prediction in the e-regime has the shape of 
a parabola with a minimum at t — T/2. This is in qualitative contrast to 
the cosh behaviour, which is standard in large volumes. £ affects both, 
the curvature and the minimum of this parabola, whereas F n only occurs 
in the additive constant — that feature is helpful for its evaluation. 

A first comparison of this curve to lattice data was presented in Ref. 
[65], using at (3 = 6, am q = 0.01 on lattice volumes 10 3 x 24 and 12 4 . 
The first among these volumes — with a linear size of L ~ 0.93 fm - 
turned out to be too small: the data for (A^(t) A 4 (0))i )2 were practically 
flat in t and incompatible with the parabola of eq. (J4.8)) for any positive 
E. This observation was consistent with the lower bound for L that we 
also found for the agreement of the microscopic spectrum with chiral 



Another observation in that study was that the corresponding history 
in v = is plagued by strong spikes, giving rise to large statistical errors. 
A huge statistics (O(10 4 ) topologically neutral configurations) would be 
required for conclusive results (see also Ref. jEE])- These spikes occur 
for the configurations with a tiny (non-zero) Dirac eigenvalue \\ t p, and 
it agrees again with chiral RMT that such configurations are most fre- 
quent in the topologically neutral sector. As a remedy to this problem, 
a method called "low mode averaging" was designed [6Tj . 

However, without applying that method we obtained a decent agree- 
ment with the prediction ()4.8|) in our second volume mentioned above 
(V ~ (1.12 fm) 4 ) in the sector \u\ = 1, and still a reasonable shape - 



RMT. 
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although somewhat flat — at \u\ = 2 [55] . In view of the leading LECs, 
it seems unfortunately impossible to extract a value of £ from such data, 
since the theoretical curvature depends on it only in an extremely weak 
way (for instance, even an extreme change from E = to £ = (250 MeV) 3 
had such a small effect that it is practically hopeless to resolve it from 
lattice data). 11 

On the other hand, F n can be extracted quite well from the vicinity 
of the minimum at t = T/2, but the value found in Ref. [HS] was too 
high. 12 

Next a study of that kind appeared in Ref. jSH], which also used D^, 
at (3 — 5.85 and fi = 1.6, now on a 10 3 x 20 lattice. These authors 
analysed the sectors \u\ = and 1 (without "low mode averaging") and 
arrived at F n = (98.3 ± 8.3) MeV. As a reason for the limitation to 
\u\ < 1, Ref. [HE] refers to the condition \v\ <C (^ 2 )- As we mentioned in 
Subsection 4.2, one expects (z/ 2 ) oc V (up to lattice artifacts), hence this 
limitation was imposed by the volume. 

Here we present again results at /3 — 5.85 on a 12 3 x 24 lattice, where 
also the latter condition admits \v\ =2, c.f. Table EH We evaluated for 
both, -DovHF and Dn, the axial-vector correlators at am q = 0.001, 0.003 
and 0.005, which turns out to be safely in the e-regime. Our propagator 
statistics at these quark masses is given in Table |U We then fitted the 
data to eq. (|4.8|) by using the chirally extrapolated factors Za (1.17 for 
-DovHF and 1.45 for D-^), along with the £ values that we obtained from 
the microscopic Dirac spectra (see Subsection 4.2 and Table 0J). For each 
of the overlap operators we performed at each of the quark masses a 
global fit over the topological sectors that we considered, which is shown 
in Figs. Uni and [T7| In particular our results for -D ov hf reveal for the first 
time a quite clear distinction between the sectors \v\ = 1 and \v\ = 2 
- this property could not be observed for D]<s up to now. For at 
am q = 0.005 we also include the neutral sector; as expected it has by far 
larger error bars than the charged sectors, but it is helpful nevertheless 
to reduce the error on F v in the global fit. 

The emerging values for the pion decay constant are consistent, 

Fir = (110 ±6) MeV (using L> ovHF ) , 

F w = (109 ±4) MeV (using D^) . (4.10) 

As an experiment, we also considered a simple re- weighting of the 

11 Only in the sector v = the sensitivity to S is significant, but there we run into 
the statistical problem mentioned before. 

12 In Ref. [HSj we gave a value around 130 MeV, but the analysis did not handle 
the renormalisation constant Za very carefully — being more precise in this aspect 
reduces the value to about 120 MeV. 
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Dirac operator 


-DovHF 


£n 


am q = 0.001 






# of propagators at \v\ = 1 


50 


50 


# of propagators at \v\ =2 


50 


50 




(110 ±8) MeV 


(109 ± 11) MeV 


am, = 0.003 






# of propagators at z/ = 1 


50 


50 


# of propagators at \v\ =2 


50 


50 


F v 


(113 ±7) MeV 


(110 ± 11) MeV 


am q = 0.005 






# of propagators at v = 




100 


# of propagators at = 1 


50 


100 


# of propagators at \v\ =2 


50 


100 




(115 ±6) MeV 


(111 ±4) MeV 



Table 4: Our results in the e-regime for the pion decay constant F WI 
based on the axial-current correlation function. The results are obtained 
at /3 = 5.85 on a 12 3 x 24 lattice. We give our statistics of the propagators 
in the different topological sectors at various bare quark masses in the e- 
regime. (Of course, different configurations where used for -D ov hf and for 
D^.) The results for were determined from fits to the quenched xPT 
formula fl^.&l ) in the range t/a 6 [11, 13], see Fias. XTb\ and\F\ 
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Figure 16: Lattice data with -D ov hf vs. predictions by quenched xPT for 
the axial-current correlation functions in the e-regime, measured sepa- 
rately in the topological sectors \v\ — 1 and 2. The global fit at each mass 
corresponds to the values of given in Table^ 
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Figure 17: Lattice data with -Dn vs. predictions by quenched xPT for the 
axial-current correlation functions in the e-regime, measured separately 
in the topological sectors \v\ — ; 1 and 2. The global fit at each mass 
corresponds to the values of given in Table^ 
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axial-current correlators involved, by means of the factor 



ni 



m„ 



A i.p 



2/i 



(4.11; 



which is (to a very good approximation) part of the fermion determinant. 
Since we take the statistics inside fixed sectors, the factor mjj"' does not 
matter here, but the second factor attaches weights to the contributions, 
which differ in particular for very small m q . As an example, we show 
in Fig. the result obtained in this way for the overlap-HF at am q = 
0.001. Of course, this is a modest step towards a 1 flavour re-weighting, 
which works well in some cases if a few hundred low lying eigenvalues are 
involved [GSJinH]- Still, in the present case we observe for the overlap-HF 
data an improved agreement with the predicted curves for \u\ = 1 and 2 
at t values relatively far from T/2. 



A 
O 

< 
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v 




Figure 18: Lattice data with -D ovH f vs. predictions by quenched xPT for 
the axial-current correlation functions, as in Fia.\lb\ but here the data are 
re-weighted with the first non-zero Dirac eigenvalue, i.e. with the factor 
\4..11\j . This re-weighting — which is most powerful at minimal m q - 
improves the agreement with the theoretical prediction at relatively large 
\t - T/2\. In this case, the fit leads to F n = (112 ± 7) MeV. 



4.4 Evaluation of based on the zero-modes 

At last we still consider an alternative method to evaluate F n in the e- 
regime. This method was introduced in Ref. [70], and it involves solely 
the zero-mode contributions to the pseudoscalar correlation function. 
Hence we now work directly in the chiral limit. Let us briefly summarise 
the main idea of this approach. 
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Ref. [70: computed the chiral Lagrangian to the next-to-next-to-lead- 

(2) 

ing order in quenched x?T, C ' PT . It can be written in a form that 
involves an auxiliary scalar field <£>o , which is coupled to the quasi Nambu- 
Goldstone field U by a new LEC denoted as K. The auxiliary field also 
contributes 

£(2)[$o] = 2X W ^° + S| $ ° (412) 

to Cr}p T , which brings in a and m as another two quenching specific 
LECs, in addition to K. The inclusion of $o supplements the quenching 
effects; in the dynamical case it decouples form the Nambu-Goldstone 
field. 

It is somewhat ambiguous how to count these additional terms in 
the quenched e-expansion. Ref. [ZO] assumes the action terms with the 
coefficients chq and K^fW~ c to be of 0(1), whereas the one with mo is in 
0(e). In particular the last assumption is a bit unusual; for instance, it 
disagrees with the framework referred to in Subsection 4.3. However, it is 
an acceptable possibility, which simplifies this approach since it removes 
the auxiliary mass term from the dominant order. If one further defines 
the dimensionless parameter 

AN 2 KF 

a = a , (4.13) 

then only the LECs F n and a occur in this order. 

For Nf valence quark flavours, we now consider the correlation func- 
tion of the pseudoscalar density P(x), which can be decomposed into a 
connected plus a disconnected part, 

V 2 (P(x)P(y)) = N^fay) - N*P 2 (x,y) (4.14) 
Pi{x,y) = Ti[i'j 5 (D + m q y 1 ( y x,y)-i'j 5 (D + m q y 1 (y,x)] 
P 2 (x,y) = Tr[ry 5 (L> + m q )~ 1 (x,x) ■ i>y 5 (D + m q )~ l {y,y)} . 

Then one performs a spectral decomposition of the propagators and ob- 
tains the residuum in terms of the zero-modes, 

m ^o(m q V) 2 (P(x)P(0)) u = N f C§(x)+Np$(x) 

connected: C^(x) = —(vj(x)v k (x) ■ vl(0)vj(0))\ u \ 
disconnected: C$(x) = (v](x)vj(x) ■ t4(°H(°)>M • ( 4 - 15 ) 



The vectors Vj denote the (exact) zero-modes of the Ginsparg- Wilson 



operator at mass zero, Dq^Vj = 0. In the terms for C, *j the zero-modes 



are summed over. 
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Next we consider the spatial integral J d 3 x P(x)P(0). Now the above 
procedure for the correlation function leads to functions CP;, (t), i = 1,2, 
which are given explicitly in Ref. [70]. In principle, these functions could 
be measured and fitted to the predictions in order to determine F n and 
a. In practice, however, it is much better to consider instead just the 
leading term in the expansion at t = T/2, 



T^Tt C ^ {t)lt=T/2 



D®s + 0(s 3 ) 



T 

t 

2 



1,2 



(4.16) 



The slopes D)A tend to be stable over a variety of fitting ranges s G 
[— Smax, Smax], s max = a, 2a, 3a.... To be explicit, the slope functions 
[771] in a volume V = L 3 x T take the form jjj] 



D 



(i) 
M 



2\u\ 



w + 



a 



01 



(F n L) 2 t 1 ' 1 ' 2.V, 
7i 1 /7 



F2VV 



+ 



(2) 



2 

2|z/| 



24 V3 v V 2 



T 2 \ 



(4.17) 



+w\ 



1 



a 



2iV c mjv 



'71-1(11 

.2 24 V 3 
where in our case 

p 1 = 0.1314565 , 



2<^> 



7i 



r 2 \ 



-E- 

1 9 ^ D i, 



1 



12^sinh^(T|p|/2) 



(4.18) 



-0.083291 



/3i is a shape coefficient, which we computed for our anisotropic volume 
according to the prescription in Ref. |44j . 

We evaluated the LECs F n and a from fits to the linear term in eq. 
fl4.16|) . We used all the zero-modes that we identified in the topological 
sectors \v\ — 1 and 2 — the statistics is given in Table El For (u 2 ) (which 
enters the expressions for through Witten-Veneziano relations) we 
inserted the result that we measured in each case, which is also given in 
Table El For each of our lattice sizes and each type of overlap operator we 
performed a global fit over both topological sectors involved, in a fitting 
range s max . The emerging optimal values for F n and a are shown in Figs. 
ITTJ1 and 1201 an d the values at s max / a = 1 are given in Table El We see 
that the results for different lattice spacings and overlap Dirac operators 
are in good agreement, and we obtain the most stable plateau for -D ov hf- 

The value that we now obtain for F n is below the one of Section 
4.4, which used a different observable and a different e-counting rule for 
the quenched terms. In fact, the result of this Section is close to the 



33 



Dirac operator 


DovHF 








5.85 


5.85 


6 


lattice size 


12 3 x 24 


12 3 x 24 


16 3 x 32 


F« 


(80 ± 14) MeV 


(74 ±11) MeV 


(75 ± 24) MeV 


a 


-17 ±10 


-19 ± 8 


-21 ± 15 



Table 5: Our results in the e-regime for the pion decay constant F w - 
along with the quenching specific LEC a — based on the zero-mode con- 
tributions to the pseudoscalar correlation function, see Subsection 4-4- 
The joint statistics in the sectors \v\ = 1 and 2, given in Tabled con- 
tributes. We give the results at fitting range s max = a, which is most 
adequate in view of eq. \4.16\ ). 



200 p 

180 

160 



12 J x 24 p=5.85 overlap-HF 
12 3 x 24 |3=5.85 standard overlap 
1 6 3 x 32 p=6 standard overlap 



140- 

> 

" 120 

LL 

100 

80- 

60 

4 L 



I 



2 

/a 



Figure 19: The results for based on a global fit of our data to a 
quenched xPT prediction for the zero-mode contributions to the pseu- 
doscalar correlations function. Here and in Fig. [H] we show the results 
of a two parameter fit over the ranges s G [T/2 — s max ,T/2 + s max ]. 



phenomenological value (we repeat that the latter amounts to ~ 86 MeV 
if one extrapolates to the chiral limit |12])- This result, as well as the 
negative value for a, are also somewhat below the values reported in 
Ref. j7D] based on the same method. Some differences are that Ref. j7U] 
always used _D N (with various values of /x), cubic volumes, a continuum 
extrapolated value for (v 2 ) and partial fits were performed. We suspect 
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Figure 20: The results for the quenching specific LEC a, based on a 
global fit of our data to a quenched X-PT prediction for the zero-mode 
contributions to the pseudoscalar correlations function. Here and in Fig. 
[TQwe show the results of a two parameter fit over the ranges s G [T/2 — 

^maxj T/2, -\- S max ] . 



that the anisotropic shape of our volumes, T = 2L, could be the main 
source of the deviation from those results [7T] . 

5 Conclusions 

We have constructed an overlap hypercube Dirac operator -D ov hf, which 
is especially suitable at a lattice spacing of a ~ 0.123 fm. It has a strongly 
improved locality compared to the standard overlap operator _D N . This 
operator defines chiral fermions on coarser lattices than D^. 

We performed quenched simulations with T> ov hf and with D N in a 
volume V ~ (1.48 fm) 3 x (2.96 fm) at (3 = 5.85 and at (3 = 6. 

In the p-regime we applied -Dovhf and measured the meson masses m n 
and m p , the PCAC quark mass mpcAc an d the pion decay constant F n 
at bare quark masses ranging from 16.1 MeV to 161 MeV. The results 
for m n and m p are similar to the values found previously with D-^ on the 
same lattice, which confirms their validity. On the other hand, mpcAc 
turned out to be much closer to m q than in the standard overlap formu- 
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lation. This implies an axial-current renormalisation constant close to 1, 
= 1.17(2), which is favourable for the connection to perturbation the- 
ory. Regarding F n , it turned out that the data obtained in the p-regime 
can hardly be extrapolated to the chiral limit. 

We considered a large number of topological charges defined by the 
fermion indices of -D ovH f or of D N (which coincide in part), and we found 
histograms which approximate well a Gaussian. The resulting topological 
susceptibility is in good agreement with the literature. 

In the e-regime we determined a value for the chiral condensate from 
the distribution of the lowest eigenvalues. For both, -D ov hf and we 
obtained £ values around (300 MeV) 3 . 

We evaluated F n in the e-regime in two ways, from the axial-current 
correlation and from the zero-mode contributions to the correlation of the 
pseudoscalar density. These two methods handle the e-counting of the 
quenched terms differently, and they yield different values for F n . The 
axial-current method leads to F n w 110 MeV, which is consistent with 
various other quenched results in the literature. The zero-mode method 
(which might be more sensitive to the anisotropy of the volume) leads to 
a lower F n around 84 MeV, in the vicinity of the phenomenological value. 
The final result of Ref. [HZ] - using again a different method, based on 
the AI =1/2 rule, still in the e-regime — is in between. We add that 
recently further methods were proposed to evaluate F n in the e-regime, 
involving As = 1 transitions [72] and a chemical potential [73*] . 

From the current results, we conclude that the methods applied here 
work in the sense that they do have the potential to evaluate at least 
the leading LECs from lattice simulations in the e-regime. The quenched 
data match the analytical predictions qualitatively (if the volume is not 
too small) and — in the setting we considered — they lead to results in 
the magnitude of the LECs in Nature. However, the quenched results 
are ambiguous: different methods yield different values. 

For precise values and a detailed comparison to phenomenology, sim- 
ulations with dynamical quarks will be needed. In particular the e- 
regime requires then dynamical Ginsparg- Wilson fermions. For instance, 
quenched re-weighting already leads to a distribution of the microscopic 
Dirac eigenvalues and the topological charges as it is expected for one 
dynamical quark flavour in view of truly dynamical QCD simula- 
tions, first tests show that it is hard to arrange for topological transitions 
|74j . but fortunately the e-regime investigations can work even in a fixed 
topological sector. Therefore, and also in view of the lattice size, the 
e-regime is promising, if one is able to handle sufficiently small quark 
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masses in a Hybrid Monte Carlo simulation, and if ergodicity inside a 
fixed topological sector is achieved. 
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